Influence of band structure eff^ects on domain-wall resistance in diluted ferromagnetic 

semiconductors 
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Intrinsic domain-wall resistance (DWR) in (Ga,Mn)As is studied theoretically and compared to 
experimental results. The recently developed model of spin transport in diluted ferromagnetic semi- 
conductors [Van Dorpe et ai, Phys. Rev. B 72, 205322 (2005)] is employed. The model combines the 
disorder-free Landauer-Biittiker formalism with the tight-binding description of the host band struc- 
ture. The obtained results show how much the spherical 4x4 kp model [Nguyen, Shchelushkin, and 
Brataas, cond-mat/0601436 overestimates DWR in the adiabatic limit, and reveal the dependence 
of DWR on the magnetization profile and crystallographic orientation of the wall. 
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Recently, current-induced domain-wall displacement^ 
and intrinsic domain-wall resistance^' have been ob- 
served in ferromagnetic semiconductor (Ga,Mn)As. 
These results attract a great deal of attention as they 
address the question of spin dynamics in the presence of 
spatially inhomogeneous magnetization texture as well as 
open the doors for novel concepts of high-density memo- 
ries and logic devices. 

Since in (Ga,Mn)As the domain-wall (DW) width 
ttXw ~ 20 nm is much longer than the mean free path 
£ ~ 0.5 nm, it has been found^ that disorder, inher- 
ent to carrier-controlled ferromagnetic semiconductors, 
plays an essential role in determining the magnitude of 
intrinsic domain- wall resistance (DWR). R can be ex- 
pected that band structure effects, particularly the spin- 
orbit interaction, which accounts for magnetoelastic^'^ 
and magnetotranspor^ properties of these materials, also 
affect the DWR. Indeed, according to recent work of 
Nguyen, Shchelushkin, and Brataas,^ in the presence of 
the valence-band splitting into light and heavy hole sub- 
bands, DWR remains non-zero even in the absence of 
disorder and in the adiabatic limit, ttXw Xf, where 
Xf is the Fermi wavelength. This conclusion, obtained 
within the 4x4 Luttinger model, questions the generality 
of the time-honored result of Cabrera and Falicov,^ who 
showed that the transition probability of the quantum 
spin through a DW approaches 1 in the adiabatic limit. 

In this paper we implement multi-orbital tight-binding 
(TB) Hamiltonian into the Landauer-Biittiker formalism 
for spin-dependent quantum transport in the ballistic 
regime. The TB method is known to be more universal 
than kp approaches. It successfully describes small quan- 
tum dots and heterojunctions between materials with dif- 
ferent band structures. This comes at the cost of a limit 
to size of the studied system. We find, however, that 
this limitation is not prohibitive in the case we study 
here. Actually, a variant of this theory was already ap- 
plied to describe the spin polarization of electron cur- 



rent in the (Ga,Mn)As/n-GaAs Esaki-Zener diodes- as 
well as to successfully evaluate the magnitude of TMR in 
(Ga,Mn)As/GaAs/(Ga,Mn)As trilayer structurespiS 

In this work we determine the magnitude of DWR for 
various domain-wall widths, shapes, and crystallographic 
orientations. In particular, we find out to what extent the 
4x4 spherical Luttinger model, ^ known to overestimate 
the effect of the spin-orbit interaction for the relevant 
values of Fermi energies;^ overvalues the magnitude of 
DWR in the adiabatic limit. Furthermore, we reveal the 
dependence of DWR on the crystallographic orientation 
of the DW, an effect not studied so far. Finally, we com- 
pare results of our disorder-free theory to experimentally 
determined values of DWR. 

The transmission matrix employed here is determined 
in terms of the extended transfer-matrix methodi''^ with- 
in the TB approach, with the spin-orbit coupling in- 
cluded. The (Ga,Mn)As valence-band structure for small 
Mn content was shown to be similar to that of GaAs.^'^ 
Thus we use sp^s* parameters proposed for GaAs^ This 
model reproduces correctly the band structure of GaAs 
in the relevant part of the Brillouin zone. Moreover, it 
takes into account the Dresselhaus terms, essential for 
the spin-dependent transport. 

The presence of the Mn ions and, therefore, of the sp-d 
exchange interaction, is introduced to the TB Hamilto- 
nian using the virtual crystal and mean-field approxima- 
tions, taking into account the appropriate weights of Ga 
and As orbitals in the wavefunctions close to the center 
of the Brillouin zonei^*i^ For the Mn content of 5% the 
hole subband splitting at the F point obtained by the 
present method is 150 meV and 46 meV for the heavy 
and light holes respectively, which compares favorably 
with the kp values of 150 and 50 meVi^ In what follows 
we use ho as a measure of magnetization, where 3/io is 
the heavy-holes splitting, and ho/2 = Bq, where Bq is 
the splitting parameter of Ref. d. 

To study coherent transport through the DW we de- 
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fine three regions along the current direction: left and 
right leads, and a central multi-layer region of the length 
Lq ^ Xw containing the DW, where Lq, taken here 
to be 40 nm, can be identified as the phase coherence 
length of the system, Lq « L^{T). The Mn concentra- 
tion is constant throughout the three sections. In the 
central region magnetization M rotates, so that its de- 
viation 9 from the initial orientation [001], taken here as 
the z-direction, changes by small steps in the consecutive 
atomic layers, according to cos9 = tanh(a;/AH'), where 
X is the distance from the DW center in the current di- 
rection. We consider the Bloch and the Neel DWsfi^ for 
which the vector M(0) x M{9) is along the current direc- 
tion or perpendicular to it, respectively. In addition to 
180° walls we study a 90° Bloch DW, for which the mag- 
netization in the left lead is normal to the M direction 
in the right lead. Since the virtual crystal approximation 
restores translational invariance in the direction perpen- 
dicular to the current, we model the lateral extent of the 
system by Born- von Karman boundary conditions. We 
assume no strain in the entire structure. 

In our calculations, we place the Fermi level Ep at ei- 
ther 0.08 eV, 0.2 eV, or 0.36 eV below the mid-point of 
the heavy and light hole bands at the F point, thus cover- 
ing the experimentally relevant hole concentration range, 
as p = 1 X 10^0 cm-3, 3.5 x 10^" cm'^, and 1 x lO^^ cm'^, 
respectively. In order to produce a non-zero current den- 
sity, we introduce a bias V by shifting the Fermi level of 
the right lead with respect to the left one. The bias value 
^ 10^^ eV is much smaller than typical splitting in the 
valence band, so that the calculated resistances are bias- 
independent. We obtain the current density from the 
transmission coefficients by energy integration over the 
range from Ep ~ eV to Ep and summation over fc|| vec- 
tors over a two-dimensional gridfi& We find that typically 
grids of 200 fc-points are sufhcient, except for more sub- 
tle problems, such as that of adiabatic passage (Fig.|21), 
which required a 550 fc-point grid. Using the current 
density and the bias we obtain the total resistance Rt of 
the structure containing the DW, in units ilfirn^. We de- 
fine DWR as a difference between Rt and the resistance 
for the uniform magnetization in the whole structure. 
Since disorder is neglected, the latter corresponds to the 
Sharvin resistance Rs, as all accessible channels are fully 
transmitting for spatially homogeneous magnetization. 

Figure n depicts DWR obtained within our sp^s* TB 
model for ho = 50 meV. We consider current flowing in 
either [110] or [100] directions, while M in the left (right) 
lead is parallel (anti-parallel) to the [001] direction. We 
present relative resistance of the DW, which is defined 
as Rr = {Rt - Rs) /Rs- In Fig. [H the DW width Xw is 
given in units of the Fermi wavelength Af of the heavy 
holes for fc^? in the [100] direction in the absence of mag- 
netization (JiQ = 0). This corresponds to Xp = 5.7 and 
3.7 nm for p = 10^" and 3.5 x 10^° cm~'^, respectively. 
The value Avi^ = denotes a step-like DW profile. 

We see in Fig. ^a) that in the adiabatic limit, Xw > 
Af, the relative resistance Rr saturates and remains non- 
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FIG. 1: (a) The predicted relative domain-wall resistance Rr 
for (Ga,Mn)As as a function of the DW width parameter Aw 
(in units of the Fermi wavelength Af ~ 5.7 nm) for various 
current directions and domain-wall profiles, (b) Rr for two 
hole densities p = 1 x 10^° cm"^ and p = 3.5 x 10^° cm"^. 



zero. The range Xw > Af is relevant for experiments 
with typical (Ga,Mn)As films. 

Some residual variations of Rr in this region are prob- 
ably caused by Fabry-Perot-like interference that oc- 
curs for a finite Lq. At the same time, according to 
the results depicted in Fig. ^b), the magnitude of Rr 
becomes significantly smaller for higher values of the 
hole density. Furthermore, Rr for the Neel wall is 
seen to be substantially lower than for the Bloch wall 
(Fig. ^a)). These findings confirm results obtained by 
Nguyen, Shchelushkin, and Brataas^ though our more 
complete theory leads to Rr up to a factor of two smaller 
than that resulting from the spherical 4x4 model. This 
can be seen for example in Fig. ^b), where the lower 
curve corresponds to ho/Ep — 0.23. In this case our TB- 
based model predicts Rr «0.14, while the kp one yields 
almost 0.3. Moreover, we find that the magnitude of Rr 
varies significantly with the DW orientation (Fig.^a)). 
This effect results from the dependence of Rt on the crys- 
tallographic directions of the current and magnetization. 

To elucidate the origin of the non-zero DWR in the adi- 
abatic limit, we have repeated our calculations switching 
off the spin-orbit coupling. Thus, we take entire complex- 
ity of the valence band into account except for putting 
the spin-orbit parameters A° and A"^ of the TB Hamilto- 
nian equal to zero in the whole structure. As shown in 
Fig. 121 Rr no longer saturates but instead vanishes with 
Xw/Xp, a behavior first noted by Cabrera and Falicov,^ 
for the free electron case. In the case of [110] direction 
the Rr exhibits certain deviation from the expected ex- 
ponential drop. We ascribe this behavior to the presence 
of grazing incidence channels for which the k component 
along the current direction is small, so that the reflection 
coefficient is still relatively high.* We conclude that non- 
zero DWR of (Ga,Mn)As in the adiabatic limit originates 
from the spin-orbit coupling rather than from multiband 
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transport in the valence band. This substantiates the 
previous suggestions^ put forward by noting that DWR 
vanishes if 4 x 4 Luttinger Hamiltonian is assumed to be 
diagonal. In this case and in the thick-wall limit, the 
spin of the carrier can adiabatically track the local mag- 
netization and the charge transport is not affected by the 
presence of DWiS*ii 
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FIG. 2; The theoretical predictions of intrinsic domain- wall 
resistance in the absence of the spin-orbit coupling. In the 
case of the domain wall normal to [10 0], the relative domain- 
wall resistance (Rr) drops exponentially to zero as a function 
of the domain- wall width parameter normalized by the Fermi 
wavelength. In the case of wall normal to [110], the grazing 
incidence channels cause a slower decay of Rr- Spin-orbit 
interaction is set to zero in this calculation, but all other 
parameters remain the same as in the calculations leading 
to a large value of Rr in Fig. The value of Af is the same 
as in Fig. Q 

In order to gain further understanding of the physics 
underlying the DW resistance, we study the thick DW 
limit. We first consider the Bloch wall, for which mag- 
netization is always perpendicular to the current, chosen 
here to flow along the [100] direction. Owing to the sym- 
metry of the zinc-blende structure, there is no anisotropic 
magnetoresistance (AMR) in such a case. The lack of 
AMR does not imply that the distribution of open chan- 
nels in the reciprocal space is identical for all magnetiza- 
tion directions 9, where is the angle of M with respect 
to the [001] direction. We show now that the depen- 
dence of this distribution on 9 gives rise to DWR in the 
adiabatic limit. To this end we use the adiabatic trans- 
port model, defined as follows. We consider all k vectors 
for which carrier propagation for given spatially uniform 
M{9) and Fermi level is possible. We project these vec- 
tors on a plane perpendicular to the current direction, 
obtaining a 2D closed area S{6) in the reciprocal space. 
Next, we turn the magnetization by a small angle and 
form the projection again, keeping constant Hq, Ep, and 
the tunneling direction. For each step, the correspond- 
ing S{6) is different, reflecting the rotation of M, and 
the coupling between the spin and orbital degrees of free- 
dom. Finally, after rotating the magnetization through 
the half-circle, we find a common subset of all S{9) for 



< 6* < TT, by retaining only these channels that are open 
for all 0. This approach is a generalization of a qualita- 
tive model presented in Ref. 0. However, we do not 
use the spherical approximation for the bands, so that 
an S{9) for a general 9 cannot be obtained by a simple 
rotation of S{0). 

In Fig. 131 we present the total DW resistance (in the 
units of Rg) as a function of 0, which is determined from 
the common subset Sc{d) of areas S corresponding to all 
relevant subbands and magnetization angles up to 9, and 
5(0) yields the Sharvin resistance Rs- We find a very 
good agreement between the calculated resistance of the 
thick Bloch DW and the one obtained in the adiabatic 
model. A direct inspection of transmission coefficients re- 
veals that the distribution of open channels is very similar 
in both approaches. To verify our arguments we consider 
a 90° DW. We find good agreement between the calcu- 
lated value of DWR and the one from the adiabatic model 
for 6* = 90°, see Fig.|3 
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FIG. 3: Total resistance normalized by the Sharvin resis- 
tance, Rt/Rs, calculated within the adiabatic model as a 
function of the angle between the magnetization and [001] di- 
rection for the Bloch and Neel walls. Horizontal lines show 
Rt /Rs computed within the full model in the adiabatic limit 
(Avk/Af ~ 1.4) for three types of walls. 

It is known that (Ga,Mn)As films reveal resistance 
anisotropy when magnetization is tilted towards the cur- 
rent directioni^ii^ The above reasoning, however, remains 
unchanged for the case of the Neel wall, and the magni- 
tudes of DWR can be obtained in the same way. In agree- 
ment with the results discussed above (see Fig.QJa)), we 
find that more channels are left open after rotating M 
by TT in the xz plane (Neel DW) than in the zy plane 
(Bloch wad). 

The values of DWR calculated from the adiabatic and 
the exact model for the 180° Neel wall differ somewhat 
more than in the case of the Bloch DW, see Fig.O This 
difference is caused by grazing incidence channels, for 
which the energy-integrated transmission coefficients are, 
on average, much smaller than the ones in the case of the 
Bloch DW. 

In the limit of vanishing spin-orbit splitting, shapes 
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and orientations of the areas S are identical for all val- 
ues. Therefore, their common subset is equal to S{0) 
and, hence, to the Sharvin resistance Rs- Thus DWR 
vanishes totally in the adiabatic approximation. How- 
ever, when the DW thickness decreases, a gradual clos- 
ing of the tunneling channels takes place. Accordingly, 
the magnitude of Rr calculated within the exact model 
grows exponentially, when Aw/Af decreases (Fig.|21l. 



1e-01 



S ie-02 

a 



1e-03 



1e-04 



Bloch Wall, y|| [ 1 1 0] 



o 



20 -3 

■ — ■ p = 1.0 x10 cm 

O— Op = 3.5x1 0^°cm"^ 
* ■ -^p^ 1.0x1 O^^cm'^ 



Mn content (%) 



FIG. 4: Predicted values of the domain-wall resistance in 
the adiabatic limit for various hole concentrations. The data 
are for the saturation value of spin-splitting at a given Mn 
concentration x. 

An important question arises whether the adiabatic 
DWR discussed here could explain experimental mag- 
nitudes of DWR in (Ga,Mn)As. Figure 2| shows the 
values of the intrinsic DWR, i.e., Rt — Rs in the adi- 
abatic limit for experimentally relevant ranges of the 



hole concentrations and Mn densities. The computed 
values of DWR are positive. The sign of the effect is 
in disagreement with the findings of Tang et al.^ who 
reported a negative value of DWR for their samples of 
(Ga,Mn)As. By contrast, Chiba et al.^ found a positive 
value, DWR= 0.5 ± O.I ilim^ for p « 2 x 10^" cm^^ and 
X = 0.05. A comparison of this experimental value with 
the results of Fig. ^ shows that the computed values are 
at least one order of magnitude smaller. We recall at this 
point that high hole and Mn concentrations, necessary to 
observe the ferromagnetism, result in a short mean free 
path in these systems, typically much smaller than the 
DW width. 

In conclusion, according to the detailed quantitative 
studies presented here, the spin-orbit interaction leads 
to a non-zero value of domain-wall resistance even in 
the adiabatic limit. The magnitude of this adiabatic 
resistance depends significantly on the character of the 
domain-wall and its crystallographic orientation. How- 
ever, the sign and magnitude of the computed results 
are in variance with the recent experimental results for 
(Ga,Mn)As. We take this disagreement as an evidence 
for the paramount importance of disorder effects in the 
physics of DWR in carrier-controlled ferromagnetic semi- 
conductors. However, non-zero domain- wall resistance 
induced by spin-orbit interactions can presumably be de- 
tected in those ferromagnetic systems, in which both co- 
herence length and mean free path are longer than the 
domain-wall width. 
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